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O ! Abstract 
D ■ 

pLn \ We give a review of the existence of Taub-NUT/bolt solutions in Einstein Gauss-Bonnet gravity 

with the parameter a in six dimensions. Although the spacetime with base space S*^ x S*^ has 
curvature singularity at r = A^, which does not admit NUT solutions, we may proceed with the 
Q-i' same computations as in the CP^ case. The investigation of thermodynamics of NUT/Bolt solutions 



in six dimensions is carried out. We compute the finite action, mass, entropy, and temperature of 



. the black hole. Then the validity of the first law of thermodynamics is demonstrated. It is shown 

^0 , that in NUT solutions all thermodynamic quantities for both base spaces are related to each other 

' by substituting = [{k + l)//c]a'^^^'^^^""^fc. So no further information is given by investigating 

NUT solution in the 5^ x S'^ case. This relation is not true for bolt solutions. A generalization 

OO , of the thermodynamics of black holes to arbitrary even dimensions is made using a new method 

O 

based on the Gibbs-Duhem relation and Gibbs free energy for NUT solutions. According to this 



^ ■ method, the finite action in Einstein Gauss-Bonnet is obtained by considering the generalized finite 

action in Einstein gravity with an additional term as a function of a. Stability analysis is done 
by investigating the heat capacity and entropy in the allowed range of a, A and N . For NUT 
solutions in d dimensions, there exist a stable phase at a narrow range of a. In six-dimensional 
Bolt solutions, metric is completely stable for B = x S"^, and is completely unstable for B — ''^^^ 
case. 
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I. INTRODUCTION 



According to the general relativity principle, the most general classical theory of gravita- 
tion in d dimensions is Lovelock gravity. It is a higher dimensional generalization of Einstein 
gravity which can extend gravity to other consistent theories such as string theory. An inter- 
esting case which is second order Lovelock gravity is Gauss-Bonnet (GB) theory with d > 5 
dimensions. In string theory context, when assuming that the tension of a string is large as 
compared to the energy scale of other variables, GB term is the first curvature correction 
to gravitation [1, 2]. Exactly the same as Einstein gravity, the Einstein-GB (EGB) field 
equation contains up to second order of metric tensor and is free of ghost. This theory has 
been extended to AdS and dS spacetimes [3]. Taub-NUT/Bolt-AdS (TN/BAdS) is a class 
of interesting spacetimes with negative cosmological constant which were investigated in 
EGB gravity [4]. The original four-dimensional solution of Taub-NUT/Bolt spacetime [5, 6] 
is only locally asymptotically flat. The spacetime develops a boundary at infinity by twisting 

bundle over S^, instead of simply being x S'^. The boundary properties of TNAdS 
spacetimes has been discussed in Einstein gravity [7]. In general, the Killing vector which 
corresponds to the coordinate that parameterizes the fibre can have a zero-dimensional 
fixed point set (called a NUT solution) or a two-dimensional fixed point set (referred to as 
a 'Bolt' solution). Generalizations to higher dimensions and higher derivative gravity have 
been done [4, 8-14]. These solutions of Einstein gravity play a central role in the construc- 
tion of M-theory configurations. Indeed, the 4-dimensional TNAdS solution provided the 
first test for the AdS/CFT correspondence in spacetimes that are only locally asymptoti- 
cally AdS [15, 16], and the Taub-NUT metric is central to the supergravity reahzation of 
the Dd-hrane of type IIA string theory [17]. It is therefore natural to suppose that the 
generalization of these solutions to the case of EGB gravity, which is the low energy limit of 
supergravity, might provide us with a window into some interesting new corners of M-theory 
moduli space. 

The AdS/CFT correspondence, is a famous conjecture that relates a, {d— l)-dimensional 
conformal field theory to d-dimensional gravity theory with a negative cosmological constant 
[18]. In light of this conjecture, similar to quantum field theory, there are some countert- 
erms which remove infinities in the gravity action. The countertcrm contains some terms 
which are functions of the curvature invariants of the induced metric on the boundary. For 
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EGB gravity, similar to Einstein gravity [19-21], in computing the action and total mass, 
one encounters infrared divergences associated with the infinite volume of the spacetime. 
Recently, the counterterm for EGB gravity was obtained for d < 9 dimensions and applied 
to static and rotating black objects [22]. Also a counterterm for charged black holes (bhs) 
in GB gravity was proposed [23]. However, an alternative regularization prescription with 
any odd-even dimensions of AdS asymptotic spacetimes for any Lovelock theory has been 
proposed [24]. This approach is known as Kounterterms and uses boundary terms with 
explicit dependence on the extrinsic curvature Kab- 

The thermodynamics of black holes in Lovelock gravity attracted a lot of attention in the 
last decade [25]. In dealing with this subject, a regularized action can be obtained by using 
the counterterm renormalization method. For Taub-NUT/Bolt black holes the area law of 
entropy is broken and one may obtain the entropy by using Gibbs-Duhem relation [26]. 
Thermodynamics of Taub-NUT/Bolt black holes is investigated in Einstein and Einstein- 
Maxwell gravity in [9, 12, 27]. By extending the dimension of spacetimes higher than seven, 
one encounters so many terms in the counterterm that the computation of thermodynamic 
quantities would be very difficult. Also a counterterm for d > 9 in EGB gravity has not 
been worked out. 

In this paper we give a new method to obtain the finite action based on the Gibbs- 
Duhem relation and Gibbs free energy for NUT solutions. According to this method, all 
thermodynamic quantities in arbitrary even dimensions can be obtained. The validity of 
the first law of thermodynamics of black holes has been tested for Taub-NUT spacetimes in 
Einstein gravity [9, 12]. We also check the vahdity of this law in EGB gravity. 

The outline of our paper is as follows. We give a brief review of the general formalism 
of EGB gravity, the counterterm method and Taub-NUT/Bolt spacetime in Section II. In 
Section III, two possible Taub-NUT/Bolt solutions of EGB gravity in six dimensions are 
reviewed. Then, in sections IV and V, all thermodynamic quantities for NUT and Bolt 
solutions are obtained and thermodynamical stability is investigated, finaly, in Section VI, 
a generalization of the thermodynamics to 2k + 2-dimensional Taub-NUT- AdS spacetimes 
in EGB gravity is made. We finish the paper with some concluding remarks. 
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II. GENERAL FORMALISM 



The gravitational action of EGB gravity in d dimension is 

^G = ^ I d''x^[-2A + R + aLoB] (1) 

where A is the cosmological constant; R, R^upai ^-nd R^i, are the Ricci scalar, the Riemann 
and the Ricci tensors of the metric g^i,] and a is the GB coefficient with dimension of 
(length)^. Lgb is the GB term, 

Lqb — R/xiyjsR'^^'^^ — ^RfxvR^^ + R? ■ (2) 

In order to have a well-defined variational principle, we must consider two surface terms, 
the Gibbons-Hawking term and its counterpart for GB gravity, which are 

4"^^ = / d'-'x^K, (3) 

AGB) _ a /■ 



= -J^ / d'^^^^iJ - ^GijK^'l (4) 

JdM 

where is the induced metric on the boundary, K is the trace of the extrinsic curvature 
of the boundary, Gij is the Einstein tensor on the induced metric, and J is the trace of the 
tensor, 

J a, = \ {2KK^Kl + K^K'^K,, - 2K,,K'^Kd, - K^K,,). (5) 

We restrict ourselves to the case of a > 0, which is consistent with heterotic string theory 
[28]. Since this term is a topological invariant in four dimensions, we therefore apply this 
term to higher dimensions which consist of only second order derivatives of the metric that 
produces second order field equations. Varying the action with respect to the metric tensor 
Qi^i,, the vacuum field equation is 

Rf,^ - ^Rgi^iy + Aqi^^ + aH^„ = 0, (6) 

where 

H^iu = 2RR^i, — AR^xR^v ~ ^R'^'Rupva + '^RJ!'^^Rvpa\ — ■^giJ.i'LGB- (7) 

We want to consider asymptotic AdS spacetime which has negative scalar curvature at 
infinity. This implies the following asymptotic expression for the Riemann tensor: 

^^^^ £2 ' 



where Ic is the effective radius of AdS spacetime in EGB gravity [22] . This parameter is 



where 



with A = — (d — — 2)/2£^. It is seen from Eq. (10) that the parameter a must have an 
upper bound as 

(c^-l)(d-2) 

In the computation of the total action Jc + + I^^ and conserved charges of EGB 
solutions, one encounters infrared divergences, associated with the infinite volume of the 
spacetime manifold. This quantities are computed by using a counterterm method which 
was proposed by Balasubramanian and Kraus [20] for AdS- Einstein gravity. The counterterm 
that regularizes the action for d < 8 solutions, for EGB gravity [22] is 



where i?, R"'^ and L^q"^ are the curvature, the Ricci tensor and the GB term (2) associated 
with the induced metric 7. Also, Q{x) is the step-function with 6 (x) = 1 provided that 
X > 0, and zero otherwise. It can be seen that as a — > {U — > 1), this action reduces to the 
familiar counterterm expression in the Einstein gravity [20]. 

Having the total finite action, one can construct a divergence-free stress-energy tensor 
by varying the total action with respect to the boundary metric 7^6, (Brown and York's 
definition of energy-momentum tensor [29]). The total finite stress-energy tensor for 0? < 8 
[22] is 

STlTab = i^G + If + lb"" + Ift'''') = Kab ' lahK + 2a{Qab ' -Qlab) 

~ ^a.(^) + ^41^(2 -U)( Ra, - 1^^,r) + ilQ {d - 6) (13) 



L ' 3 ' d-3 ' ' V 2 



i T', , „ , . -\ ^nh( Rr/jR'^ — ^, . \-R'^) — ^. . \-RRnh ~l~ 2R'^R, 

^{d-3y{d-5)^ 2'""^ 4{d-2) ' 2{d-2) 



cadb 



2(d-2) " " 2(^-2)'"" ' 2(d-3)(d-4)(d-5) 
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where [30] 

Qab = 2KKacK'^ - 2KacK'"'Kab + K^biK^aK'''' -K^) + 2KRab + RKab - 2K''''R,adb - ^RacK, 

(14) 

and Hab is given by (7) in terms of the boundary metric ^ab- 

Prom the finite stress-energy tensor, one can compute conserved charges. To do this job, 
we must choose a spacehke surface E in the boundary dM. with metric Uij and write the 
boundary metric in the Arnowitt-Deser-Misner form, 

-iabdx^'dx^ = -W^di^ + (Tij{d^p' + V'dt)(dif' + V^dt), (15) 

where A/" and V* arc the lapse function and shift vector, respectively. The coordinates 
i = 1, . . . ,d — 2 are the angular variables, parameterizing the hypersurface of constant r. 
The conserved charge associated with a killing vector is 

0(0 = ^ d''-\^u''Tabi\ (16) 

where a is the determinant of the metric and u"- is the normal to quasilocal boundary hy- 
persurface E. For example the conserved mass M, is the charge associated with the timehke 
killing vector ^ = dt- 

The Taub-NUT/Bolt-AdS metric is one of the solutions of EGB gravity which was in- 
vestigated in {2k + 2) dimensions [4] . This metric is constructed on a base space endowed 
with an Einstein-Kahler metric Sg. The Euclidean section of the {2k + 2)-dimensional 
Taub-NUT/Bolt spacetime can be written as: 

ds^ = F{r){dT + NAf + F-\r)dr^ + (r^ - N^)'E.b (17) 

where r is the coordinate on the fibre and A is the Kahler form of the base space B, 
which is proportional to some covariantly constant two-form. Here N is the NUT charge and 
F{r) is a function of r which is obtained by solving the EGB field equation (6). The solution 
will describe a 'NUT' if the fixed point set of the U(l) isometry d/dr (i.e. the points where 
F{r) — 0) is less than 2/c-dimension and a 'Bolt' if the fixed point set is 2A;-dimensions. 

III. SIX-DIMENSIONAL SOLUTIONS 

In this section we give a revision of the six-dimensional Taub-NUT/bolt solutions(17) of 
GB gravity will be done. The function F{r) for all nonextremal choices of the base space B 
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can be written [4] in the form of 
where 

= 1 + pzriv^)! 

+ - + ISATV + 57V«) (19) 

and the function P{r) depends on the choice of the curved base space B. For the case 
B = CP^ the Kahler form A and the CP^ metric are 

^2 = 6sin^6(#2 + sin2 6#i) (20) 

dSa^ = 6{d6^ + sin^ ^ cos^ 6(#2 + sin^ ^i^i)^ 

^sir?i2{di^^ + sin^ cos' ixdi)^^)}. (21) 

The function P(r) in this case is 

16a2(r4 + 6r2Ar2 + AT^) 
^CP^ (^2 _ ;v2)4 ■ l^^) 

The other possibihty is B = -S' x 5"', where -S' is the 2-sphere with = d^^ + sin' 6*^0', 
and the one-form A is 

yl = 2cos^i(i0i + 2cos^2C?02- (23) 

In this case P(r) is 

32a2(r^ + 4r'7V' + 7V4) 
^^'x^' (^2 _ jv2)4 • (24) 

For asymptotically AdS spacetimes [A < provided that |A| < 5/(12Q;inax)], the function 
F{r) is not real for all values of r in the range < r < oo and a > 0. 
The Ricci scalar and Lqb in the bulk (1) are 

1 d 



R 



^^(r' - + 4rF(r)(r' - N^) - ^(r' - SA^^) 



(25) 



(r2 - Ar2)2 dr 

Las = [3F(r)(r' + AT') - (r' - AT')] + 3rF(r)' - 2rF(r) + Fr} , 

(26) 

where F = 2/3 for ^ = CP' and "1" for B = x respectively 



7 



A. NUT Solutions :ln brief, the conditions for having a NUT solution, which have 
been described completely in some previous works [9, 12], are as follows: 

(i) F{r = N) =0; (ii) F'{r = N) = l/(3iV); (iii) F{r) should have no positive roots at 
r > N. The first condition comes from the fact that all extra dimensions should collapse 
to zero at the fixed point set of d/dt. The second condition allows one to avoid a conical 
singularity with a smoothly closed fiber at r = A^. Using the first condition, the GB gravity 
may admit NUT solutions with CP^ and S'^ x S"^ base spaces when the mass parameter m 
is fixed to be 

1 fi 

^(CF) ^ _i^^(3Ajv^ + 57v2 _ 5a), (27) 

mf^ ^ --^N{6AN'^ + 10N'^ -15a). (28) 

By inserting m„ from Eqs. (27) and (28) in (18), the function -F„(r) can be obtained. But as 
it has been mentioned in Ref . [4] , the metric with B — x has a curvature singularity at 
r — N. Although the second condition to have a NUT solution is satisfied numerically, by 
having a curvature singularity at r — N, this condition is basically violated. Therefore, the 
metric does not admit a NUT solution in this case. In Fig. 1, the function F„(r) is plotted 
versus r for a = 1/35 and A = —10 for both base spaces. As it is shown in this figure, at 
the small vicinity of r = A^, -F„(r) is complex with a continues real part. Following section 
IV, we will show that if one does not consider this fact and investigate the thermodynamics 
of NUT solutions for spacetimes with both base spaces, no obvious difference will be seen. 
In the hmit of a — > 0, Fn{r) for both base factors goes to F^ut in Einstein gravity [9]. The 
allowed range of F{r) in the CP^ case is more extended than S"^ x S'^ case. Also if a goes 
to amax, or -u ^ 0, Fn{r) would be complex in some range of r so that in both base spaces, 
we find that Fn{r) is not real at a = amax for all values of r in the range of < r < oo and 
a > 0. This fact is true for NUT solutions of arbitrary even dimensions. 

The boundary of this kind of spacetime will be an asymptotic surface at some large 
radius r. The boundary metric 7,^ diverges at infinity. By rescaling upon a conformal 
factor ^c/r^ (see Eq.(9)), the dual boundary metric hab — lim£^/r^7ab converges and the 

r— »oo 

general line element of the dual field theory for even-dimensional TNAdS spacetimes on the 
boundary can be obtained as 

dsl = {dr + NA)^ + il^B. (29) 
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FIG. 1: Fn{r) versus r for B = (dotted line) and B = x ("bold-line" for pure real and 
"solid- line" for the real part in the vicinity oi N = 1). 

In fact, lim£^/r^F/v[/T(r) = 1. As it has been shown in (29), the Gauss-Bonnet coefficient 
a enters the metric on the boundary due to £c- 

B. Bolt Solutions: The conditions for having Bolt Solutions are (i) F{r — rb) — 
and (ii) F'{rb) — 1/{3N) with Vf, > N. These Conditions give the following relation for mi, 

-A _ ^ 15iVV2[4 _ A(^2 _ 3^2)] _ 40^,(^2 ^ ^2) 



rriL = 



3A{rf + 5N^) - 10(n 



I5n 



(30) 



and an equation for with the base space CP^ 

SNAn^ + (2 + 3AN^)rb^ - N{4: + 3AN^)n - 3AN^ - 6N^ + 8a = (31) 

which at least has one real solution. This real solution for N < ^/a yields rb > N. 
For the case of B = S'^ x S"^, the mass parameter is 

(5) _ 3A(rg + 5A^^) - 10(4 - 3A^^) + 15A^Vg[4 - A(rg - 3N^)] - 60a(rg + N^) 



ml 



15n 



(32) 



and rb can be found by solving the following equation: 

3A^Arb^ - 6N(AN^ + l)r6^ + 2(r^ - N^)rb + S(AN^ + 2)N^ + Aa{2rb - 3N) = 0. (33) 

This equation will have two real roots \i N < N^ax- At N — N^ax there will be only one 
rb > N. Therefore in order to have Bolt solution in this case, there is a maximum value for 
N [4]. It is found that the allowed range of r^, for this case is more extended than in the 
CP^ case. 
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IV. THERMODYNAMICS OF 6-D NUT SOLUTIONS 



A. CP^ case: Prom relations (25) and (26), we can simply compute the bulk action (1). 
By computing the other actions for the boundary and counterterm (3), (4), and (12), the 
finite action /fin = + -^6^^ + ^b^^^ ~ ^ct'^^^ can be obtained as 

« = -^(2AiV' + 5iV"-10c). (34) 

In Eq. (44), j3 = An/ F'[N) = 127rA^ is the inverse of the Hawking temperature Th which 
is found by demanding regularity of the Euclidean manifold (17). The total mass can be 
calculated as 

AlAr = 27rm„. (35) 

The total angular momentum is zero, which is the same as Einstein gravity. The Gibbs free 
energy and entropy may be computed as 

^^^^ " /3 " 2332807r4r| + 9Th' ^ ^ 

4"^) ^ _(|^) ^ -327r27V2(2AiV^ + 3N' - 2a), (37) 

Also the entropy can be obtained, as the Einstein case [9] , by using Gibbs-Duhem relation, 
which is 

S = f3(M - iiA) - I (38) 

where /Xj is chemical potential and Cj is the conserved charge (such as mass, angular mo- 
mentum,...). The Smarr type formula for mass versus extensive quantity S may be written 
as 

M^S^'\S) = -^^(3AZ^ + 5Z- 5a) (39) 
where Z = N'^ is a function of S, which may be obtained from the following equation: 

S + 327r^Z(2AZ2 + 3Z- 2a) = 0. (40) 

The Hawking temperature Th may be calculated as 

(41) 

ds 12nN ^ ' 

Now the specific heat can be obtained as 

C},"') = Th{^) = 1287r2Ar2(3AAr4 + 2>N^ _ «). (42) 
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Note that the entropy become negative for a. < 3/2 A^^ — |A| A^^ and specific heat is negative 
for a > 3{N^ - |A| N^). Thus, for thermally stable solution for B = CP^ case, the value of 
a parameter must be in the range of 

3/2Ar2 - |A| AT^ < a < 3(A^2 - |A| A^^), ^43^ 

provided that a < a^ax, N > and A < 0. In the limit of a — > 0, (Einstein gravity), this 
spacetime is completely unstable [12]. 

B. S^xS'^ case: As in the previous case, the finite action, Gibbs free energy, entropy, 
and specific heat may be obtained as 

CT' = -^(2A^^ + 57V^-15a), (44) 

^^^^ " ~ 2332807r4r5 + 3Th ^ ' 

gis^xs^) ^ -S27r^N^{2AN'^ + 3N^ - 3a). (46) 
^ifxs') = i927r2A^2^2AA^^ + 2A^2 _ (47) 

By comparing all quantities in the B — CP^ and the S'^ x 5"^ case, we can see that by 
substituting a'^^'^ = {3/2)a^'^^^'^ in all equations in the previous case, all thermodynamic 
quantities can be obtained in this base space. By expanding computations in 2A; + 2 dimen- 
sions, it can be seen that 

k 

According to this fact, there is no difference between these two base spaces and if we do 
not consider the metric with 5"^ x 5"^, as a NUT solution of EGB gravity, nothing would be 
missed. 

It is worthwhile to mention that, for all thermodynamic quantities which are calculated 
in these cases, the first law of thermodynamics, dM. — TndS, is satisfied. 
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V. THERMODYNAMICS OF 6-D BOLT SOLUTIONS 



A. S'^xS^ case: The finite action can be calculated as 



.(52x52) _ 47r^ 



{(15AiV^ + 30A^^ 



+30 NW - 70 A^Vb^ - 10 iV^A Tft^ - 8 A^^A + 10 r^^ 
+3 n^K) + (120 A^Vft^A - 300 A^^ + 100 r^^ - 120 NW 
-200A^VA- 120AiV^ + 72Ar6^)a + 480(Ar2+ rft^ja^}. (49) 

In the limit of a — > the action (49) goes to the corresponding action in Einstein gravity 
[9, 12]. The entropy and specific heat may be written as 

~ 3[2(r2 - N^) - A{rl - N^f + 4a) ^^^^^ 

+24 n^rl + 9 - 46 r^) + 6(3 + r^) - 42 nV^ + 18 nV^ 

+96 {n^ + r-Da"^ + (-84 - 24 A - 168 nV^ - 24 A 

+60r^ + 88n2Ar^ - 168n^Ar^)Q;}, (50) 

8Tr\rl -N^ + Aa) 



(5^x52) 



3[Q;(8Ar2 - 8AAr2 _ 4) + ^(r^ - mY][A{rl - N^)^ - 2{rl - - Aa] 
{ - 192 + (576 n^K N^-32A + 336 A^^ + 96 A A^^ + 144 r^) 

+ (20 rb^A^ - 12 A^A^^ + 376 A^Vb^A + 48 AT^^r^^A^ - 120 A A^^ 

-144 A^Vfe^A^ + 192 + 96 r^N"^ - 72 A^Vt^A - 120 n^k - 192 A^^ 

+88 A^V6^A2)q; - 88 A^^A r^^ + 104 A^^A + 36 A^^ + 36 + 36 A^^A 

+74rb^A=^Ar^ - 30 A^iVVfe^ - 15 A^A^Vb^ - A'ir^^t^N'^ - mr^N^ 

-36 n^N'^ - 28 Tb^A + 5 n^^k^ + 9 A^A^^o - 24 r^^A N^} (51) 

In Fig. 2 the entropy and specific heat are plotted for A — —10 and a = 1/30 as a 
function of A^ in the allowed range of A^ < A^max = 0.12, for the outer horizon of the Bolt 
solution. As it can be see from this figure, the Bolt solution in this case is completely stable. 

B. CP^ case: The finite action, entropy, and heat capacity can be obtained as 

^ 5[3A(r2 - N^y - 6(r2 - AT^) - 8a] ^ ^ ' 4Uiv j a 

-40 A^Vb - 32 A^VbA + 15 A AT^ + 45 A^^rt^A + 30 A^^ 

-10 - 15 + 60 n^N'^ + 3 Tfe^A} (52) 
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FIG. 2: Entropy (dot) and specific heat (bold-line) versus of Bolt solution for B = S"^ x S'^. 



5[8a(3Ar2 - 3Am - 2) + 3A(r2 - N^y][3A{rl - N^f - Q{rl - N^) - 8a] 
{ - 180 n^A + 450 N^n^K - 630 N^n^K + 180 rjJ'N - 300 iV^A 

+420 N'^K + 270 A^V^A - 90 K^N^n^ - 120 n^K - 360 N^n^ + 

+180 + 180 + 180 A^^A + 45 A^A^^^ + 312 N^Vb^A^ + 180 N^Vb 

-264 A^Vfe^A^ - 90 A^A^Vfe^ - 171 n^A^A^^ + 18 n^^A"^ + 72 A^^r^A^ 

+288r6^A2Ar^ - 180 rb^N^ - ISOn'^N^ - 90 rb^N A - UO^'^N^A'^ 

+10 (128 A^^ + 48A Tft^ + 96A^r6 + 128A^V6A + 96r6^A A^^ + 48A A^^) 

+ (-90 A^A^^ + 480 NWA - 640 N^n - 840 A^^ + 54 n^A^ - 480 N^b^A'^ 

-360 Nvb^A + 1520 rb^A^^A + 960 Vb^N + 720 Vb^N^ + 360 NW^^ 

-600 A AT^ - 504 N\b^A^ + 180 A^V^^A^ + 840 A^V^^A - 96 A^V^A^ 

-336 Tb^A + 576 A^V^^A^ + 760 - 904 N''rbA)a - 640 a^}. (54) 



Fig. 3 shows the entropy and specific heat for A —10 and a = 1/30 as a function of N 
in the allowed range. It is shown that the Bolt solution in this case is completely unstable. 

VI. THERMODYNAMICS OF (2A; + 2)-DIMENSIONAL NUT SOLUTIONS 



,(CP2) 




+15 A A^^ + 45 A^Vb^A + 30 A^^ - 10 r^^ - 15 A^V^^A 
+60 n^N^ + 32 iVVfoA + 3 n^A + 40 iVVj 



(53) 



,(CP2) 



167r2(r^-Ar2 + 4a) 
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FIG. 3: Entropy (dot) and specific heat (bold-line) versus N of Bolt solution for B 



In this section a new procedure to achieve thermodynamic quantities by using the Gibbs- 
Duhem relation and Gibbs free energy will be introduced. 

The generalized form of even-dimensional TNAdS spacetimes in EGB gravity has been 
constructed [4]. The mass parameter in d — 2k+2 dimensions for B — CP*^ can be calculated 
as 



r(2A; + 3) 



-(r(^))^ 



[2 A N\k + lf + {2k'^ + 3k + l)kN^ + (1 - 4 k'^) {k - l)ka\ . 



(55) 



As it has been mentioned in 6D NUT solutions, all quantities go to the Einstein gravity case 
in the limit of a — > 0. Also the GB parameter a, has not been appear in the calculation 
of the total mass and inverse Hawking temperature /5. The similar case has been shown in 
Ref [31]. Therefore we may consider the total mass and /3, similar to d-dimensional NUT 
solution in Einstein gravity [12] which have been obtained as 



= 4(A; + l)nN, 



fc-1 



m 



(56) 
(57) 



What about finite action? Up to now, the main method for computing the finite action has 
been the counterterm method, which has been used in Sec. IV for six dimensions. In order 
to extend this method to higher dimensions, one needs more terms than (12). However 
recently, a new counterterm action has been achieved for 7 < d < 9 in EGB gravity which 
has so many terms [22] . According to this method, the existence of the counterterm action 
is necessary for computing the thermodynamic quantities. Now we give a different method 
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based on the Gibbs-Duhem relation (38) and Gibbs free energy for NUT solutions that can 
calculate regularized action independent of the counterterm method. 

According to this method, first, first write the finite action in 2A; + 2 -dimensional Einstein 

gravity [12] with an additional term as a function of a. This term must be considered in 
such a way that in the limit of a 0, it goes to zero for consistency to Einstein gravity. 
Therefore the finite action in EGB gravity can be written as 

[{2k^ + 3k + 1)N^ +{k + l)2AN^ + e{k- 2) fk{a, N)] . (58) 

After that, calculate the Gibbs free energy as a function of Th and then calculate the entropy 
Sg similar to Eq. (37). The other way of computing the entropy, which is independent of 
the previous method, is to use the Gibbs-Duhem relation (38). By calculating the entropy 
in the latter way {Sod), the first order differential equation Sqd — Sq gives us fk{o(, N) as 

h{a,N) = -k{Ak^-l)a. (59) 

The results in 6D may be considered as an initial condition. The function (59) is independent 
of N and is the same as the last term in (55). Now the finite action may be calculated in 
EGB gravity by substituting (59) in the action (58). For k = 2 (QD), the action (58) is equal 
to (34) in section IV. In the appendix, computations for 8D are done using the counterterm 
method and we see that there is a good coincidence with our prescription. The Gibbs free 
energy, entropy, and heat capacity may be obtained as 

G%{Th) = ^('' + '^) r (-1/2 - k) , , (5e + 2e + Ak + l) 
' ((A;+l)Tj,)''+V*^+5/24'=+3^ ^ ^ 

-128 {4k^ + 12k'^ + llk^ + k^ -3k- 1) n^Tn^ka + A], (60) 

= -4*^-V'=-i/2)^2fe-2p(^/2-A;)r(A; + l) 

[(2 k^ + k- 1)N^ + {k + 1)2 A A^^ - (4A;^ - 8A; + 3)ka] , (61) 

= 4(fe-i/2)^fe-i/2^2fc-2p (^/2 - A;) r (A; + 1) 

[(2 k'^ + k- l)kN^ + (1 + 2A;)2 A A^^ - (4 A;^ - 12 A;^ + 11 A; - 3 )ka\ . (62) 

The spacetime for B — CP*^ is stable in the following range: 

Ar2[2(A; + l)Ar2A + (2A;2 + A; - 1)] A^2[2(A; + i)2iv2A + A;(2A;2 + A; - 1)] 

A;(4A;2 - 8A; + 3) A;(4A;3 - 12A;2 + llA; - 3) ' ^ ' 
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provided that < a < ccmax (H), N > 0, and A < 0. According to all calculations in this 
section, the validity of the first law of thermodynamics of black holes can be simply verified. 



VII. CONCLUSION 



In this paper, a review of the existence of non extremal asymptotically AdS Taub- 
NUT/Bolt solutions in EGB gravity with curved base spaces has been discussed. The 
spacetime with base space S'^ x S'^ has a curvature singularity at r — N. This fact violates 
the second condition for having a NUT solution. Therefore, the metric basically does not 
admit NUT solutions. But, as this condition has been satisfied numerically, we may proceed 
with computations the same as in the CP^ case. We show that in the limit of a ^ 0, F„(r) 
for both base spaces, goes to Fnut in Einstein gravity. Also if a goes to amax, or -u — 0, 
Fn{r) would be complex in some range of r so that in both base spaces, wc find that Fn^r) is 
not real at a — amax for all values of r in the range < r < oo and a > 0. This fact is true 
for NUT solutions of arbitrary even dimensions. The boundary of Taub-NUT-AdS space- 
time at some large radius r is obtained, and it was shown that the Gauss-Bonnet coefficient 
a enters the metric on the boundary due to £c. The investigation of the thermodynamics of 
NUT/bolt solutions in six dimensions is carried out. By obtaining the regularized action and 
stress tensor, the total mass is derived. The Smarr-type formula for the mass as a function 
of the extensive parameter S is calculated. By computing the Hawking temperature, the va- 
hdity of the first law of thermodynamics is demonstrated. It is shown that in NUT solutions 
all thermodynamic quantities for both base spaces are related to each other by substituting 
Q,cp'= ^ _|_ iy^^s^xs^x...sl_ rpj^jg relation is not true for bolt solutions. So, no further 
information is given by investigating NUT solutions in the S*^ x S'^ case, and wc can remove 
this metric from NUT solutions. A generalization of the thermodynamics of black holes in 
arbitrary even dimensions is made using a new method based on the Gibbs- Duhem relation 
and Gibbs free energy for NUT solutions. According to this method, the finite action for 
TNAdS spacetimes in EGB gravity is obtained by considering the generahzed finite action 
for TNAdS spacetimes in Einstein gravity with an additional term as a function of a. This 
term should vanished if a — 0. The great importance of this method is the computation 
of the finite action, which is easier than the counterterm prescription, for some spacetimes 
higher than seven dimensions in higher derivative gravity. However, this is just a conjecture 



16 



that its validity must be checked as an open problem. 

Finally, we perform the stability analysis by investigating the heat capacity and entropy 
in the allowed range of a, A, and A^. For NUT solutions in d dimensions, there exists 
a stable phase at a narrow range of a for both base factors. For bolt solutions in six 
dimensions, the metric is completely stable in the allowed range of a andiV for the base 
factor B — S'^ X S'^, and is completely unstable for B — CP^ case. An interesting open 
problem would be to compare the results of different prescriptions (the counterterm method, 
Kounterterm, and our procedure) in computing the finite action. 



APPENDIX: FINITE ACTION FOR EIGHT-DIMENSIONAL TNADS 
SPACETIMES IN EGB GRAVITY 

The Ricci scalar and Lgb ior d — 8 in the bulk are: 



R 



d 



(^2 _ ]y2y dr 
12 



ar ^ ^ ^ ^ 5 



Lgb 



d 



(^2 _ j^2y dr 



dF (r) 
dr 



dF (r) 



r2\2 



F (r) (5r" -2r^N^ -3N^)- (r^ - N^) 



+2 rF (r)' (5 - N^) - A rF (r) (r' - AT^) + 1 r ( - 3A^') 



In This case counterterm has been obtained [22] as 



'ct 



Stt 



— / d''-^x^{ - 



dM 

e^Q {d - 6) 



d-2 2 + U. icOid-A) 



2{d-3)^{d-5) 



U(RabR' 



4 

ab 



■)( 



3 ' 

d-1 
A{d-2 



2{d-3) 



2 - U)R 



-R' 



d-3 



(U - l)Lgg 



+0 {d - 8) 



where 



1 + 31/30(C/- 1) ]Le- 



( 3d-l 



57600 



2{d-A) 

19£^ I- \ 

^''^ {U - 1)£[-) 



j^T^abr. {d-l) (ci+1) 3 

(d-3)3(d-5)(ci-7) V4(d-2) 16(rf-2)2 

J o 

— ZK K riacbd+Tm V — ^ V rLab + 



2{d-2) 



2{d-2y 



and the third Lovelock term is 



(64) 



(65) 



(66) 



RV^R], (67) 



£(*^) o ryabcd p pe/ , o pab pee pd/ i O/l TDCibcd p pe 
(3) — ^-n- -n-cdeZ-K a6 + °-"' ed-"- 6/-"- ae + ^^ii Kcdbe^a 

+ ZRR"'^'^'^ Rcdab + '^^R"^'^'^ RcaRdh + 16i?"^i?6c-Ra ~ V2RR"'^ Rg}, + i?^ . (68) 
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All quantities in (66), (67) and (68) must be calculated in the induced metric ^ah- By 
considering all actions in the Bulk, boundary, and counterterm, the finite action can be 
obtained as 

167r2Ar3/? 



Cn) = 33 (8AiV^^ + 28iV2 - 105a). (69) 
It has a good coincidence with the action (58) for /c = 3 in Section (VI). 
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